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ABSTRACT 



Aims. Measurements of the Sunyaev-Zel'dovich (hereafter SZ) effect distortion of the cosmic microwave background provide us 
with an independent method to derive the gas temperature of galaxy clusters. In merging galaxy clusters the gas distribution is 
inhomogeneous and, therefore, the method of temperature measuring based on the SZ effect should be more relevant than that based 
on an X-ray emission analysis. Here we study a method for measuring the gas temperature in merging clusters by means of the SZ 
effect. 

Methods. Our calculations of intensity maps of the SZ effect include relativistic corrections considered within the framework of the 
Wright formalism and utilize a cosmological numerical simulation of a merging galaxy cluster evolved with its baryon physics. 
Results. We found that the gas temperature in merging clusters can be measured by means of the ratio of the SZ intensity at a 
low frequency (128 GHz) to that at a high frequency (369 GHz). This SZ intensity ratio permits us to reveal prominent features of 
the temperature structure caused by violent merger shock waves. Therefore, measurements of the ratio of the SZ intensities are a 
promising tool for measuring gas temperature in merging galaxy clusters. 
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1. Introduction 

Galaxy clusters are large structures in the Universe, with radii 
of the order of a megaparsec. The space between galaxies in the 
clusters is filled with low-density (10 _1 — 10 3 cm 4 ) high tem- 
perature (kftT ~ 1 — 10 keV) electron-proton plasma (for a re- 
view, see e.g., Sarazin 1986). Inverse Compton scattering of hot 
free electrons in clusters of galaxies on the cosmic microwave 
background (CMB) radiation field causes a change in the in- 
tensity of the CMB radiation towards clusters of galaxies (the 
Sunyaev-Zel'dovich effect, hereinafter the SZ effect; for a re- 
view, see Sunyaev & Zel'dovich 1980). 

The SZ effect is important for cosmology and the study of 
clusters of galaxies (for a review, see Birkinshaw 1999). It mea- 
sures the pressure of an electron population integrated along 
the line of sight as long as free electrons are non-relativistic. 
Relativistic effects are significant for high temperature plasmas 
in galaxy clusters (see, e.g. Rephaeli 1995). A relativistically 
correct formalism for the SZ effect based on the probability dis- 
tribution of the photon frequency shift after scattering was given 
by Wright (1979) to describe the Comptonization process of soft 
photons by mildly relativistic plasma. 

Relativistic effects for the SZ effect permit us to measure 
the temperature of intracluster plasma (see Pointecouteau et al. 
1998; Hansen et al. 2002). This method is more promising for 
measuring the temperature of merging and/or distant clusters of 
galaxies than those based on studies of an X-ray spectrum anal- 
ysis (see Pointecouteau et al. 1998). This is because the X-ray 
emission traces the denser component (the X-ray emission is 
proportional to J nl(l)dl (where « e is the electron number den- 
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sity and dl is the integrated line of sight), while the SZ effect is 
proportional to J n e (l)dl) and the cluster's X-ray surface bright- 
ness strongly decreases with redshift while the SZ brightness is 
independent of redshift. 

Using hydrodynamic simulations of galaxy clusters, Kay et 
al. (2008) compare the temperatures derived from the X-ray 
spectroscopy and the SZ effect. Their method for measuring the 
SZ temperature is based on an equation (see Eq. 6 from Kay et al. 
2008) which is only valid in the Rayleigh- Jeans limit. As known 
relativistic effects on the CMB intensity distortion are more 
significant at higher frequencies, Colafrancesco & Marchegiani 
(2010) conclude that to obtain detailed information about the 
cluster temperature distribution one must use high-frequency 
spectral observations of the SZ effect in the range 300-400 GHz 
and show that the SZ temperature can be extracted even for cool 
non-merging galaxy clusters, such as the Perseus and Abell 2199 
clusters. 

In this paper, using the relativistically correct Wright formal- 
ism and results of a cosmological numerical simulation of a cool 
distant merging cluster, we show how the SZ temperature can 
be derived by means of multi-frequency observations of the SZ 
effect. We calculate intensity maps of the SZ effect at different 
frequencies and propose to use the ratio of the SZ intensities at 
two frequencies to derive the SZ temperature of galaxy clusters. 

The layout of the paper is as follows. We describe the cosmo- 
logical numerical simulation of the cool distant merging cluster 
in Sect. 2. We calculate the SZ intensity maps at different fre- 
quencies in the framework of the Wright formalism in Sect. 3. 
We consider the ratio of the SZ intensities at two frequencies to 
find a convenient method for observing the SZ temperature of 
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galaxy clusters in Sect. 4 and present our discussions and con- 
clusions in Sects. 5 and 6. 



2. Numerical simulations 

The galaxy cluster simulations of Dubois et al. (2010) which 
we use in this paper are run with the Adaptive Mesh Refinement 
(AMR) code RAMSES (Teyssier 2002). The evolution of the gas 
is followed using a second-order unsplit Goduno v scheme for the 
Euler equations. The Riemann solver used to compute the flux at 
a cell interface is the acoustic solver using a first-order MinMod 
total variation diminishing scheme to reconstruct the interpo- 
lated variables from their cell-centered values. Collisionless par- 
ticles (dark matter, stars and black hole particles) are evolved 
using a particle-mesh solver with cloud-in-cell interpolation. 

The simulations are performed using a re-simulation (zoom) 
technique: the coarse region is a 128 3 grid with M DM = 2.9x 10 10 
M DM resolution in a 80 h^J Mpc simulation box, where /170 
is the Hubble constant in units of 70 km s Mpc~' . This region 
contains a smaller 256 3 equivalent grid in a sphere of radius 20 
h~Q Mpc with M DM = 3.6 x 10 9 M Q DM resolution, which in turn 
encloses the final high resolution sphere with radius 6 h^j Mpc, 
512 3 equivalent grid and M DM = 4.5 x 10 8 M Q DM resolution. 
The maximum level of refinement reached in this simulation al- 
lows to resolve a minimum spatial scale of 1 . 19 hrlkpc. 

A flat ACDM cosmology was assumed with total matter den- 
sity Q.,„ = 0.3, baryon density Qb - 0.045, dark energy den- 
sity Qa = 0.7, fluctuation amplitude at 8 h~ ' Mpc <x 8 = 0.90 
and Hubble constant Ho=70 km s _1 Mpc 1 that corresponds to 
the Wilkinson Microwave Anisotropies Probe 1 year best-fitting 
cosmology (Spergel et al. 2003). 

The resimulated region tracks the formation of a galaxy clus- 
ter with a 1:1 major merger occurring at z = 0.8. This z = 0.8 
major galaxy merger drives the cluster gas to temperatures twice 
the virial temperature thanks to violent shock waves. The pro- 
jected mass-weighted temperature map is shown in Fig.Q] where 
x-direction is orthogonal to the merger plane. 

Standard recipes for atomic cooling, heating from UV back- 
ground and star formation in high density regions are included. 
This simulation also includes a novel recipe for AGN feedback 
based on a kinetic -jet energy input that prevents the cooling 
catastrophe to occurs, the formation of unrealistic massive cool- 
ing flows, and regulates the star formation in the central galaxy 
(see Dubois et al. 2010 for further details on the simulation). 



3. The SZ effect from the simulated cluster 

In this section, the SZ intensity maps are calculated in the 
framework of the relativistic Wright formalism for the simulated 
galaxy cluster. 

The CMB intensity change produced by the SZ effect by 
non-relativistic electrons considered within the framework of 
the Kompaneets approximation is (see for a review, Birkinshaw 
1999): 



AI nr (x) = Iog(x)y gas 



(1) 



where x = hv/kbT cm \>, Iq = 2(k\,T cm b) 3 /(he) 2 , and the spectral 
function g(x) is given by 



x exp(x) / exp(x) + 1 

\x 4 



(exp(x) - l) 2 \ exp(x) - 1 



(2) 




Fig. 1. The mass-weighted temperature map (in keV) of the sim- 
ulated cluster along the x-direction at z=0.74 



The subscript 'nr 1 denotes the fact that Eq. was obtained in 
the non-relativistic limit. The Comptonization parameter y gas is 
given by 



m,c 



(3) 



where the line-of-sight integral extends from the last scattering 
surface of the CMB radiation to the observer at redshift z=0, 
T e is the electron temperature, « gas is the number density of the 
gas, cr T is the Thomson cross-section, m e the electron mass, c 
the speed of light, k\, the Boltzmann constant and h the Planck 
constant. 

The CMB intensity change in the Wright formalism can be 
written in the form proposed by Prokhorov et al. (2010), and is 



AI(x) = h 



0~T f 
2 H Z 



s kbT e G(x, T e )dl, 



(4) 



where the spectral function g(x) is changed to the generalized 
spectral function G(x, T e ) which depends explicitly on the elec- 
tron temperature. 

The relativistic spectral function G(x, T e ) derived in the 
framework of the Wright formalism is given by 



x 3 exp(-3s) 



exp(xexp(-s)) - 1 exp(x) - 1 



\ds(5) 



where ®(r e ) = k\,T e /m e c 2 , and Pi(s,T e ) is the distribution of 
frequency shifts for single scattering (Wright 1979; Birkinshaw 
1999). 

There are three basic spectral features that characterize the 
thermal, non-relativistic SZ effect signal: a minimum of its in- 
tensity is located at a dimensionless frequency xi = 2.26 (v=128 
GHz), the crossover frequency is xo = 3.83 (v=217 GHz), and a 
maximum of its intensity is located at a dimensionless frequency 
X2 = 6.51 (v = 369 GHz). 

Observations of the SZ effect close to the crossover fre- 
quency are biased due to the presence of the kinematic al SZ 
effect (for a review, see Birkinshaw 1999), which is associ- 
ated with the peculiar velocity of the galaxy cluster. The un- 
known value of the peculiar velocity limits the ability to measure 
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the cluster temperature directly through the displacement of the 
crossover frequency of the SZ effect in the correct relativistic 
treatment (Colafrancesco et al. 2009). 

Prokhorov et al. (2010) show that the choice of frequencies 
x\ = 2.26 and xn = 6.51 corresponding to minimum and max- 
imum values of the SZ intensity in the Kompaneets approxi- 
mation is suitable to analyze mildly relativistic electron popu- 
lations. Since the frequency of its maximum is located at 369 
GHz, observations at this high frequency should be promising to 
derive the temperature by means of the SZ effect (Colafrancesco 
et al. 2010). 

To produce the SZ intensity maps at frequencies x\ = 2.26 
and xn — 6.51 we use the 3D density and temperature maps for 
the simulated galaxy cluster considered in Sect. 2. We calculated 
the SZ effect using the Wright formalism. The intensity maps of 
the SZ effect at these frequencies derived from the simulation 
maps of the gas density and temperature are plotted in Figs. [2] 
and[3] respectively. 




Fig. 2. The intensity map I/Iq of the SZ effect at a frequency 128 
GHz derived from the numerical simulation in the framework of 
the Wright formalism. 



The morphologies of the SZ intensity simulated maps at fre- 
quencies X] = 2.26 and X2 = 6.51 are similar. This is because 
the SZ effect from cool galaxy clusters can be approximately 
described in the framework of the Kompaneets approximation. 
Note that the spectral function g(x) does not depend on gas tem- 
perature. However, we now show that the gas temperature struc- 
ture of a cool galaxy cluster can be derived by means of a more 
accurate analysis. 



4. The ratio of the SZ intensity at frequency 128 
GHz to that at frequency 369 GHz 

In this section, we show that the ratio of the SZ intensities at 
frequencies 128 GHz to 369 GHz provides us with a convenient 
method for measuring the SZ temperature of galaxy clusters. 

Multi-frequency analysis of the SZ effect permits us to de- 
rive temperature maps for galaxy clusters. Colafrancesco & 
Marchegiani (2010) fitted the six frequency SZ effect data (v = 
300, 320, 340, 360, 380 and 400 GHz) with a relativistic model 
of the thermal SZ effect and assuming, for each experimental 




Fig. 3. The intensity map I/Iq of the SZ effect at a frequency 369 
GHz derived from the numerical simulation in the framework of 
the Wright formalism. 

data point, an uncertainty of 0. 1 %. They present results of the fit- 
ting procedure used to extract the cluster temperature from a set 
of simulated spatially resolved spectroscopic SZ effect observa- 
tions in different bands of the spectrum. The Perseus and Abell 
2199 clusters studied in Colafrancesco & Marchegiani (2010) 
are relaxed galaxy clusters, here we present our results of an 
analysis of the temperature structure of the simulated merging 
galaxy cluster. 

Using the relativistic Wright formalism for the SZ effect we 
found that the ratio of the relativistic spectral function G{x\, T e ) 
at frequency 128 GHz to that of function G(x2, T e ) at frequency 
369 GHz is a monotonic function of temperature. For an isother- 
mal galaxy cluster the ratio of the SZ intensities is given by 
AI{xi)/AI(x 2 ) = G(x u T e )/G(x 2 , T e ). Note that the ratio of the 
SZ intensities does not depend on temperature in the framework 
of the Kompaneets formalism. 

We have checked the obtained monotonical dependance us- 
ing the generalized Kompaneets equation derived by Challinor 
& Lasenby (1998) including relativistic effects. Using their Eq. 
(28) which is valid for k^T e < 10 keV we find that the ratio of 
the SZ intensity at frequency 128 GHz to that at frequency 369 
GHz equals 



AZ(*i) 
A/(x 2 ) 



-0.607- 2.143 x 



?n e c 2 



(6) 



Comparing Eq. © from this paper and Eq. (6) from Kay et 
al. (2008) we find that the temperature derived from the ratio of 
the SZ intensities at frequency 128 GHz to that at frequency 369 
GHz within the framework of the Wright formalism is equiva- 
lent the Compton-averaged electron temperature for cool galaxy 
clusters for which Eq. (28) from Challinor & Lasenby (1998) is 
valid. 

The ratio of the SZ intensity maps at frequencies 128 GHz 
and 369 GHz derived within the framework of the Wright for- 
malism (see Figs. [2] and [3} is shown in Fig. [4] We note that 
regions of high temperature on the mass-weighted temperature 
map (see Fig. Q]) corresponds to regions with low values of the 
SZ intensity ratio as expected from Eq. ©. 

Using the monotonical dependance of the SZ intensity ratio 
on temperature derived within the framework of the Wright for- 
malism we derive the SZ temperature map from the map of the 
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Fig. 4. The ratio of the SZ intensity at frequency 128 GHz to that 
at frequency 369 GHz for the simulated cluster 



ratio of the SZ intensities at frequencies 128 GHz to 369 GHz. 
The SZ temperature map is shown in Fig. [5] 




Fig. 5. Temperature (in keV) of the simulated cluster along the x 
direction at z=0.74 derived from the ratio of the SZ intensities. 



The temperature maps in Figs. Q] and [5] show high temper- 
ature "arc-like" structures. The origin of these "arc-like" struc- 
tures on the temperature maps is the major merger of two galaxy 
clusters occurring at z=0.8. The evolution of the mass-weighted 
temperature along the x-direction during its merging phase is 
shown in Fig. 5 from Dubois et al. (2010) at different redshifts 
z=0.88, z=0.74, and z=0.58. As shown by Dubois et al. (2010) 
this major merger drives the cluster gas to temperatures twice 
the virial temperature thanks to violent merger shock waves. 

Figure [5] shows that the SZ temperature is higher than the 
mass-weighted temperature (see Fig. [TJ. This is consistent with 
the results of Kay et al. (2008). 

We find, from Eq. (6), that a sensitivity of order of ~ 0.3% of 
the predicted SZ effect is necessary to measure the SZ tempera- 



ture in cool galaxy clusters with a precision of 1 keV by means 
of the SZ intensity ratio. 

Our analysis shows that the ratio of the SZ intensities is a 
promising approach for measuring the temperature structure in 
merging galaxy clusters. 

5. Discussions 

A source of bias in the observations of the SZ effect could be 
provided by a possibly relevant kinematic SZ effect. The kine- 
matic SZ effect arises from the bulk motion of the medium rela- 
tive to the CMB rest frame. In this section we propose a method 
to extract the kinematical SZ effect from SZ observations and 
also discuss various constraints on the ability of measuring the 
gas temperature in distant merging galaxy clusters by means of 
X-ray observations. 

5. 1 . The kinematic SZ effect 

To study the kinematical SZ effect from SZ observations, 
Rephaeli & Lafav (1991) proposed to measure the SZ effect at 
the crossover frequency. The contribution of the kinematical SZ 
effect is maximal at the crossover while the thermal SZ effect 
in the Kompaneets approximation equals zero. Using Eq. (28) 
from Challinor & Lasenby (1998), we calculate the ratio of the 
intensities of the relativistic SZ effect to the kinematic SZ effect 
at the crossover frequency. This ratio is given by 

/„!(* = 3.83) ^ / k h T e \ 2 300km/s 
Iun(x = 3.83) ~ ' Ukevj v 



(7) 



where v is the peculiar velocity of a galaxy cluster and 300 km/s 
is the rms value of the peculiar velocity distribution of galaxy 
clusters (see Giovanelli et al. 1998). 

Equation|7]shows that the relativistic corrections to the ther- 
mal SZ effect in galaxy clusters at the crossover frequency can be 
significant and biases measurements of peculiar velocity by the 
Rephaeli & Lafav method, particulary in hot clusters with tem- 
perature ~ 10 keV. Therefore, below we propose another method 
to measure peculiar velocities. We note that the approximate 
function described the SZ relativistic corrections taken from 
Challinor & Lasenby (1998) for k B T e < 10 keV has two fre- 
quencies at which relativistic corrections equal zero. The values 
of these frequencies are x a = 3.33 and x\, = 8.02. Measurements 
of SZ intensities at these frequencies provides us with a method 
to derive the peculiar velocity. Using Eq. (28) from Challinor & 
Lasenby (1998), we find 



v m e c z AI{x a )g(x b ) - M{x h )g(x a ) 

n e al | — = x 

c Iqctj hix^gix^) - Mx^gixa) 



(8) 



where the spectral function h(x) describing the kinematic SZ ef- 
fect is 



h(x) = 



1 exp(x) 



(exp(x)-l) 2 - 



(9) 



The optical depth of the gas to Compton scattering can be ob- 
tained from spectral and spatial X-ray observations, and can be 
determined accurately in relaxed galaxy clusters. The proposed 
method to derive the peculiar velocity of a galaxy cluster based 
on Eq. (O is not affected by biases due to relativistic corrections 
of the thermal SZ effect. 
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It is important that measurements of the value of (J n e dljv/c 
permit us to extract contributions of the kinematical SZ effect 
from SZ observations at frequencies of 128 GHz and 369 GHz, 
since the contribution of the kinematical SZ effect is proportional 
to (^J n e dtj v/c. 

5.2. Constraints on the ability of measuring the gas 
temperature via X-ray observations 

To compare methods for measuring the gas temperature in merg- 
ing clusters based on the SZ effect with that based on X-ray 
emission, we produce the X-ray surface brightness map of ther- 
mal bremsstrahlung emission for the simulated cluster. For the 
sake of illustration, the normalized X-ray surface brightness map 
of the simulated cluster in the [2.0-10.0 keV] band in logarith- 
mic scale is shown in Fig. [6] which can be obtained by modern 
X-ray satellites, such as XMM-Newton, Chandra, and Suzaku. 
The region corresponding to the highest temperature in Fig. |5]is 
shown by a black circle in Fig. [6] The X-ray surface brightness 
of this region is two orders of magnitude smaller than the maxi- 
mal surface brightness value in Fig. [6] while the SZ intensity (see 
contours of the SZ intensity superimposed on the X-ray surface 
brightness map, the inner contour represent 10% of the maxi- 
mum of the SZ intensity and the outer contour represents 1 % of 
the maximum of the SZ intensity) of this region is one order of 
magnitude smaller than the maximal SZ intensity value in Figs. 
[2] and [3] Therefore the X-ray emission traces the denser compo- 
nent while the SZ effect provides us with a method to study more 
rarified gas regions. 
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Fig. 6. The normalized X-ray surface brightness map of the sim- 
ulated cluster in the [2.0-10.0 keV] band in logarithmic scale. 
The region corresponding to the highest temperature in Fig. [5]is 
shown by a black circle. 

The Abell 3376 cluster at redshift z=0.046 is a cool nearby 
cluster and provides us with an example of galaxy clusters 
where a shock has not been revealed by X-ray observations with 
the XMM-Newton X-ray satellite, but instead it has been de- 
tected by radio observations with the Very Large Array telescope 
(Bagchi et al. 2006). This is because the X-ray surface brightness 
in the region of the east ringlike structure is very weak and does 
not permit us to detect this shock. 



Since the X-ray surface brightness decreases with redshift as 
1/(1 + z) 4 , the X-ray surface brightness for the simulated clus- 
ter at z=0.74 decreases by a factor ~9 compared with the same 
cluster at z=0. The decrease of the X-ray surface brightness with 
redshift constrains the ability of measuring the gas temperature 
in distant galaxy clusters by means of X-ray observations. 

6. Conclusions 

Merging galaxy clusters are an interesting astrophysical labo- 
ratory for studying gasdynamic processes. Mergers of galaxy 
clusters are very energetic astrophysical events in which huge 
gravitational energy is released. In the course of a merger, a sig- 
nificant portion of this energy, which is carried by the gas, is 
dissipated by merger shock waves. This leads to a heating of the 
gas to higher temperatures. Studying the temperature structure of 
merging galaxy clusters is important to reveal heated gas regions 
associated with merger shock waves. 

Studying X-ray spectra provides us two independent ap- 
proaches to determine the gas temperature in galaxy clusters, this 
is because the bremsstrahlung continuum spectrum depends on 
the gas temperature and ratios of emission line fluxes are func- 
tions of the gas temperature (e.g. the ratio of the He-like to H- 
like iron line fluxes, see e.g. Prokhorov et al. 2009). An analy- 
sis of the SZ effect permits us to determine the gas temperature 
by studying deviations of the SZ intensity spectrum from that 
derived within the framework of the Kompaneets approxima- 
tion. Colafrancesco & Marchegiani (2010) show that the method 
based on the SZ intensity deviations can be interesting for mea- 
suring the gas temperature even in cool clusters, such as the 
Perseus and Abell 2199 clusters, if uncertainties of observational 
data are 0.1% - 1%. In this paper, using the relativistically cor- 
rect Wright formalism and results of 3-D hydrodynamic numer- 
ical simulations of a cool distant merging cluster, we show how 
the SZ temperature can be derived by means of multi-frequency 
observations of the SZ effect. 

To produce a realistic cool distant merging cluster we use a 
zoom cosmological simulation from Dubois et al. (2010). We 
find that the simulated cluster undergoes a violent merger at 
z=0.8, which drives the cluster gas to higher temperatures. The 
temperature map of the simulated cluster becomes very inhomo- 
geneous because of the merger activity (see Fig.[TJ. We consider 
a method for measuring the temperature structure in this merg- 
ing simulated cluster by means of distortions of the SZ intensity 
map due to relativistic effects. 

We incorporate the relativistic Wright formalism for model- 
ing the SZ effect in the numerical simulation using the algorithm 
proposed by Prokhorov et al. (2010). We calculate the SZ inten- 
sity maps at low (128 GHz) and high (369 GHz) frequencies and 
find these SZ intensity maps look similar as expected for a cool 
galaxy cluster. To provide a method for measuring the tempera- 
ture structure in cool merging galaxy clusters we propose to use 
the ratio of the SZ intensities at two frequencies. 

We calculate the ratio of the SZ intensity at a low frequency 
(128 GHz) to that at a high frequency (369 GHz) for the sim- 
ulated merging cluster and show that this ratio is a promising 
method to measure the SZ temperature in galaxy clusters. The 
calculated map of the SZ temperature shows that main features 
(such as "arc-like" structures) of the temperature structures of 
merging clusters may be revealed by using this method. 

The next generation SZ effect experiments discussed in 
Colafrancesco & Marchegiani (2010) are needed to reach the 
required high sensitivity with the purpose of being independent 
in measuring gas temperature in galaxy clusters. Experimental 
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configurations which are basically the same as those of the 
Millimetron project (see http://www.asc.rssi.ru/millimetron/) 
will have the power to measure the gas temperature via the SZ 
effect in any cluster, including cool ones. 
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